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Abstract 

I carry out a finite-size scaling study of the correlation length in SU(3) lattice gauge theory 
coupled to 12 fundamental flavor fermions, using recent data published by Fodor, Holland, Kuti, 
Nogradi and Schroeder [1]. I make the assumption that the system is conformal in the zero-mass, 
infinite volume limit, that scaling is violated by both nonzero fermion mass and by finite volume, 
and that the scaling function in each channel is determined self-consistently by the data. From 
several different observables I extract a common exponent for the scaling of the correlation length ^ 
with the fermion mass niq, ^ ~ ruq with ~ 1.35. Shortcomings of the analysis are discussed. 
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A recent paper by Fodor, et al [l|] presents an analysis of spectroscopy for SU{3) gauge 
theory coupled to Nf = 12 flavors of fermions. This theory is a potential candidate for 
beyond - Standard Model physics. However, it is a subject of some recent controversy. 
Appelquist et al [2I, 0|, performing a calculation of a running coupling constant, concluded 
that in the limit of vanishing fermion mass, it had an infrared-attractive fixed point (IRFP). 
In that limit it exhibits conformal behavior at long distances. The authors of Ref. [l| 
collected spectroscopic data at one value of the gauge coupling, four simulation volumes, 
and eight fermion masses, for a total of twelve volume - mass combinations. They analyzed 
their data under the competing assumptions that the system was confining and chirally 
broken, or conformal, and concluded that their data favored the confining and chirally 
broken scenario. Other references relevant to this controversy include (4-13 



The "conformal scenario" for a system like the one under discussion assumes that the 
long distance behavior of the theory is described by one relevant coupling, the fermion mass 
rriq. In infinite volume, tuning the mass to zero causes the correlation length to diverge 
algebraically, 

i_ 

^ ~ mg . (1) 

The quantity i/m is the leading relevant exponent for the system, in statistical physics lan- 
guage. This exponent is related to the anomalous dimension 7^ of the mass operator ipip, 
and determines the running of the mass parameter according to 

= -7..(/)m(/i), (2) 

Um = i + Imig*)- All other couplings, including the gauge coupling (more properly, the 
distance of the gauge coupling from its fixed point value) are irrelevant. 

However, no simulation is ever done in infinite volume. The system size L is also a 
relevant parameter since the correlation length only diverges in the 1/L — )■ limit. When 
the correlation length measured in a system of size L (call it becomes comparable to L, 
C,L saturates at L even as ruq vanishes. However, if the only large length scales in the problem 
are ^ and L, then overall factors of length can only involve ^ and L. For the correlation 
length itself, this argument says that 

a = LF{^/L) (3) 

where F{x) is some unknown function of ^/L. A somewhat more useful version of this 
relation invokes Eq. [H to say 

a = Lf{Ly-mq). (4) 

Then one can plot ^l/L vs Ly^^mq for many L's, and vary y^. Under this variation, data 
from different L's will march across the x axis at different rates. The exponent can be 
determined by tuning to collapse the data onto a single curve. 

The authors of Ref. [Ij confronted a subset of their data with Eq. [T]and all their data with 
Eq. [3l The reason for this note and my analysis is that they chose a particular functional 
form for F{x) in their fit. However, in general, the actual form of the scaling function is 
unknown. Limiting behavior is known: for example, at large x, F{x) = x and at small x, 
F{x) goes to a constant. This means that finding F{x) or f{x) is itself part of the fit. In 
addition, there is no reason for different observables to have the same F{x). The finite box 
size can - and the data show that it does - affect them differently. 
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Not having a-priori knowledge of the scahng function means that it is difficult to assign 
a goodness-of-fit parameter, like a chi-squared, to a determination of i/m- All one can do is 
to compare the |/m's from different data sets, and ask if they are consistent. It also means 
that the resulting value of ym will have a large uncertainty. However, I think it is still a 
potentially informative task, to ask, whether the data of Ref. [l| is consistent with the finite 
size scaling hypothesis, while letting the data itself determine the scaling function. That is 
the subject of this note. 

This methodology was used in Ref. to measure ym{g'^) in the SU{3) - Nf = 2 sextet 
system. It produced relatively noisy exponents. The technique of using correlation functions 
in the Schrodinger functional is much more accurate, but the finite-size scaling exponents 
agreed reasonably well with these better measurements (l5| . 



The technique has already been described in Ref. [IJ], so we will just proceed to results. 
I have analyzed the mass spectra of the following states from the data sets of Ref. [ij : the 
pseudoscalar (would-be Goldstone boson in a chirally broken theory), vector and axial vector 
mesons, baryon, and pseudoscalar decay constant, I will define the correlation length C,l 
to be just the inverse mass, (or 1/fn) in a lattice whose spatial length is L. 

To begin the analysis, we have to see if the data shows the appropriate qualitative be- 
havior: does seem to fiatten out at small ruq, at an L - dependent value? Fig. [1] shows 
that (with one exception) the trend is as expected, and furthermore larger correlates with 
larger L. Data at smaller masses on larger volumes would be desirable (was that not always 
so?) to push to the independent regime, but in addition, smaller volume data would 
do as well. The size of the finite volume effect is, not surprisingly, different for different 
observables. 

The exception is /tt, panel (b), which has a tiny L dependence with an unexpected order, 

falls with L, and no plateau yet observed, as f-,, always monotonically decreases with niq. 
The axial vector matrix element (0|y4o|vr) ~ mTtfj, does approach a plateau because does. 

By eye, before the different L data separate, the correlation lengths seems to show the 
power law behavior of Eq. [1] Nevertheless, the figure illustrates the danger of a simple fit 
to a power law: if the system is conformal in the zero mass limit, the smallest fermion mass 
data is presumably closest to conformality, and yet it is the most contaminated by finite 
volume effects. While the data does not show this, presumably the largest masses could 
be far enough away from the critical region that they might be outside it, following some 
different scaling law. One is then forced to consider cutting the data from both the high 
and low mass ends, to produce a fit to Eq. [H not a desirable procedure. 

Next we perform a scan of ^l/L vs L^""mq, varying ym, using the cleanest data set, the 
pseudoscalar. This is shown in Fig. [21 for ym = 1.0, 1.2, 1.4, 1.6. Collapse to a single scaling 
curve seems to be occurring. Notice that ^l/L is not zero; the volume is not infinite. 

Finally, I attempt to determine a best fit value of ym- The method is that of Bhattacharjee 
and Seno [l6|. The idea is to use each data set (each different L values, for the same channel) 
to estimate the scaling curve and to find the ym which pulls the other L sets onto it. This is 
done inclusively; all data sets take a turn at being the fiducial. The quantity to minimize is 

p(yj = E E E - /.(ir™..)) ' (5) 

over . / - \ J 

p jf^p I, over ^ J ' 

Data set p is used to estimate the scaling function j[x). This is done by interpolation, either 
by polynomials or rational functions, using the recorded values of ^l/L. The label "over" 
indicates that the sum only includes data from set j whose x values, L^^rriij, overlap the 
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FIG. 1: Correlation length versus inverse quark mass. Panels use the inverse pion mass (a), /^r 
(b), and proton mass (c). Plotting symbols are for different simulation volumes, squares, L = 48; 
diamonds, L = 40; crosses, L = 32; octagons, L = 24. 

range of x's of set p. The overall factor of 1/Nover counts the total number of points and 
guards against recording a zero value of P if there are no overlap. P is minimized by the 
optimal Um- This is folded into a jackknife. 

Again, the scaling function is not identical for different particle correlators, and so I 
choose not to combine the data further. All data sets tested so far produce similar values for 
Umi about 1.35. Results are given in Tabled with errors from a single-elimination jackknife. 
Bhattacharjee and Seno advocate taking an error from an approximation to the second 
derivative of P, 



This gives slightly smaller uncertainty estimates, about 0.1 at r/ = 1. The scaling curves for 
this Um are shown in Fig. El 

The authors of Ref. [ij confronted their data with Eq. [T]and concluded that the hypothesis 
was disfavored. They repeated their analysis using Eq. 121 making a specific choice for F{x), 
and reached a similar conclusion. The choice of a particular functional form for F{L/^) is, 
I have already remarked, unjustified, and the authors' analysis may simply show that their 
functional form for F{x) is not the one which actually describes the data. 

My analysis of the data of Ref. [l] with the assumption that the infinite volume theory 
is conformal while the conformality is broken by both the quark mass and the finite volume 
produces a consistent picture, that the leading relevant exponent at their simulation param- 
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FIG. 2: Plots of vs niqLy^ with S^l = l/''7^7r for four choices of ym- (a) y™, = 1.0, (b) 

= 1-2, (c) ym = 1-4 (d) y^ = 1-6. Plotting symbols are for different spatial sizes, squares, 
L = 48; diamonds, L = 40; crosses, L = 32; octagons, L = 24. 



eters is ym = 1-35 or 7^ = 0.35, with unfortunately unimpressively large uncertainty. Of 
course, this is an analysis for which the data set was not designed. It could be improved by 
more mass values at all chosen volumes. 

The small 'jm measured here resembles results from other nearly-conformal theories ob- 
served to date jisl, 17-191 . 

Notice, finally, that this is far from being a complete story. For the fit itself, one could 
be concerned with, and include, non-scaling contributions. (See Ref. j20| which does this. 
The authors specified their F{x) rather than letting the data do so.) 

More importantly, the ym which comes out is very likely not to be an actual scaling 
exponent. A few moment's reflection shows why: Because the gauge coupling runs so slowly, 
simulations done over a small range of volumes cannot flow to a fixed point (if it exists) 
unless they begin very close to it. 

That this is expected, is easy to see from the one-loop beta function result, where under 
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FIG. 3: Plots of ^l/L vs nigLy"^ for ym = 1.35: (a) pseudoscalar (would-be Goldstone) (b) /tt, 
(c) proton. Plotting symbols are for different spatial sizes, squares, L = 48; diamonds, L = 40; 
crosses, L = 32; octagons, L = 24. 



a scale change s the inverse coupling is shifted by 

1 1 h 



log s + . . . 



(7) 



In three-flavor SU{3), bi = 9 and in 12 flavor SU{3), hi = 3, so that the equivalent scale 
changes in the two theories, for an equal change in coupling, are Si2 = In ordinary QCD, 
the coupling runs from weak at a distance of 0.1 fm to strong at a distance 1.0 fm, or over 
s = 10. It is hard to get a large aspect ratio s from a set of numerical simulations at a 
single set of bare parameters. Therefore, whether or not a running coupling in one of these 
many-fermion theories actually has an IR fixed point, it runs so slowly that for all practical 
purposes its running can be neglected. Then the zero mass limit is effectively conformal. 
That is all that is needed to motivate a scaling analysis as is done here. 

(The lowest order beta function argument is merely suggestive of problems, but if the 
theory has an IRFP, than the situation on the weak coupling side calls for even slower 
running than the one loop formula.) 
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TABLE I: Exponent Um from various hadronic channels. Errors are from a single-elimination 
jackknife. 

channel ym 

pseudoscalar 1.35(23) 

nucleon 1.43(26) 

U 1.23(31) 

vector 1.33(22) 

axial vector 1.32(12) 



Taking this argument further, it says that a measurement of Um at one set of bare cou- 
phng values does not answer the question of whether the theory actually has an IRFP. 
Simulations which have made predictions for an exponent map out in some prescription 
and ^m{g^), both from simulations at many values of the bare parameters. They then sep- 
arately determine the critical coupling (in some scheme) gl and read off 'jmigl)- (Compare 



Refs. [15|, ll7H19l|.) Most likely, finite size scaling studies of spectroscopy are just not the 
method of choice for discovering whether the theory has an IRFP, and if it does, accurately 
determining either g^ or ■yrnigl)- 
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